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The equation of the original curved surface of end harmonic gearing is 
determined by displacement analysis of flat flexible gear. The dis- 
placement analysis is also used to calculate the strength and rigidity 
of the gear. The latter is regarded as a circular plate with two 
concentrated loads, since its torsional rigidity is much larger than 
its bending rigidity. Small-deflection theory of thin plates is used . 
to solve for the displacement, of any point in the middle plane of 
the gear. New expressions are given for given for radial and tangential 
displacements of the middle plane under asymmetrical loading. A digital 
computer is used to jbtain numerical values for the displacements. 
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TM 774*3 

A DEFORMATION ANALYSIS OF FLAT FLEXIBLE GEAR AND ITS 
EQUATION OF ORIGINAL CURVED SURFACES 
SHEN YUNWEN* 

ABSTRACT /175** 

Based on the elastic thin plate theory, the displacement of 
an arbitrary point on the flexible gear affected by a wave 
generator was investigated, an equation for the original curved 
surface which is required to study the engaging theory of the end 
harmonic gearing was derived, analytical solutions as well as 
calculated results of "u", "v" , and "w" were obtained. Thereby, 
this paper offers the theoretical basis for studying the 
geometric theory and the strength calculations of the end 
harmonic gearing transmission. 

FORWARD 

Although it has been twenty years since the beginning of the 
development of the harmonic gearing transmission, yet current 
research and application of the harmonic gearing are limited to 
radial engagement transmission. Along with the extended 
application of the harmonic gearing transmission, some equipment 
require more strict requirements for the axial dimensions as well 
as the torsion rigidity of this type of transmission. 

Consequently, end surface engaging of harmonic gearing 
transmission is being developed. In foreign countries, some 
patents have been awarded. However, due to the complexity of the 
theoretical study, few theoretical papers with systematical 
analyses have been published. Moreover, some papers even 
presented dubious results. For example, reference (1) presented 
an erroneous result with an imperfect analysis. Reference (2) 
used excessive approximations. Due to the urgent applications of 
this transmission method, a systematical study on the end 
harmonic gearing transmission was conducted. This paper not only 
pointed out the error in reference (1), but also presented an 
analytical method for the deformation of a flat flexible gear. 

*Associate Professor, Northeastern Polytechnical University 
** Numbers in margin indicate foreign pagination. 
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I 

In addition, it also Introduced the original curved surface 

equation which is necessary for the engaging analysis of the end 

harmonic gearing transmission. 

II. MATHEMATICAL MODEL OF DUAL WAVE TRANSMISSION 

When we study the flat flexible gear of the end harmonic 

gearing transmission, we can treat the flat gear as a thin plate 

of uniform thickness which is held along the circumference of an 

inside circle with a radius Cr = r ). 

c 



Figure 1. Diagram of the Flat Flexible Gear 

As shown in Figure 1, two concentrated loads p are applied on the 
free boundary of the thin plate which also takes a torsional 
moment "T" in the middle plane. For more precise analysis, the 
effect of the bending moment due to engaging which will be 
distributed according to a certain pattern in the engaging zone 
should also be considered. For simplification, we temporarily 
ignore the effect of this engaging force. This effect on the 
stress distribution will be solved by the proper adjustment of a 
certain constant (Ref. 3). 

From the plate and shell theory (Ref. 4), we can see that 
the stress and the deformation state can be described by the T. 
Karmen' s non-linear equations. Let us assume that the middle 
plane stress function is $ , and the displacement of any point on 
the middle plane along the z-axis is w, then in a polar 
coordinate ( r * < ?), the group of non-linear equations can be 
written as 


/ lie 


2 


dw 


*♦*£)♦(* 

a. i +3(1 *<J> \ d ( 1 \1 , 9 

r 1 d<P l ) dr 1 dr \r dr ) dr \_r dr ) J + 3 1 

♦A £)}. 


dw 

dr 


( 1 ) 


where, D - cylindrical rigidity of the plate D= E6^/12( 1-u^) ,E and 
m, elastic modulus and Poisson ratio respectively; 5 - plate 
thickness, i.e., the wall thickness of the flat flexible gear; q 
- uniformly distributed load on the plate in the lateral 
direction; V* - operator, 

= * + ± _*_ + l a ‘ 

dr 1 r dr + r 1 d<P l • 


From the engineering point of view, it is not necessary to 
include the geometric factors and mechanical factors to solve 
these equations, because these two factors would complicate the 
problem. This problem can be simplified into a small deflection 
problem because the effect of the deformation due to the torsion 
moment 'T M can be neglected since the rigidity of the flat 
flexible gear is high. In addition, the gear is in a very slight 
bending configuration because the maximum axial deformation W 
of the flexible gear under the concentratedload is about the same 
as the working height of the gear-tooth or the gear thickness. 

Since a circular plate is under a concentrated load "P", the uniform 

distributed load can be considered as q=0. When we substitute 4> 

= 0 , then equation ( 1 ) can be simplified as 


/ & Id 1 d l \fd l w l dw 1 d^iv \ _ (2) 

v v “ = (‘a?‘ + 7 sr+p v 1 wr 11. 

Let u, v, and w tepresent the radial, tangential, and axial 
displacement respectively, based on the condition that the 
strain of the middle plane is zero, and according to the 
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elastic thin plate theory we can easily find the geometric 
equations as follows 


V 


du 

dr 


i / 9w y 

+J\bt) 


0» 

du> 


dv • 1 / 1 du> \* 

‘•‘7ft* y*i\T -ft) 

; 1 du dv t dw \( 1 du; \ 

l x - w 7 -ft*sr*{sr)\T -ft)-°- 


(3) 


The physical equations are found to be [4] and [5] 


and 


i# r , /I dm . 1 

-d[- S t + «(7 -ar + r' 


(4) 


- radial, tangential, and shear strain 


(5) 


e and e m 
<p rq> 


/ 1 7 7 


M and M - the bending moment and torsion 


where, e r 
respectively; M 

moment along the directions of r and 9 respectively; and - 
shear forces along the directions of r and ? respectively. 

III. AXIAL DISPLACEMENT 


The axial displacement <*» can be solved with equation (2). 

The following sections will find the solutions of these 
differential equations. 

1. Boundary Conditions 

For a flat flexible gear with an inner holding circle r = r c 
the boundary conditions of the inner holding side can be 
expressed as 


(w) rm , e = 0, 



( 6 ) 
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Since the outside circumference (r = r^) of the flexible gear is 
free, hence, 

(A/ r ) r . r> «0i (7) 

In addition, there are two concentrated loads, 180° apart from 
the x-axis of the free boundary; therefore the shear force is not 
continuous. When the positions of the two concentrated forces 
are treated in such a way that the phase difference is x , then 
the period will be formed by the accumulation of a 2x - cycle 
function. In the meantime, if we assume that the effect of the 
concentratedload is equivalent to a uniformly distributed load q 
acting on the minute arc 2r^.A<p, then we can use the method in 
Reference (6) to expand the concentrated load into a Fourier 
series, and use relation of q^ =P/2r« A * to determine the 
corresponding constants. Finally, let us find the limiting value 
by assuming A9"*0, and accumulate the effect of the two 
concentrated forces which have a phase difference of x. Then the 
expanded Fourier series will be 

p f = -2/L(-l+ E cos nr), 

xr, V 2 » .2.4.6. ' 


and we can also find the other boundary condition of the free 
boundaries as 



rc<p t, mrt xr, \ 2 


E 


cos nr 


)• 


(8) 


2. Solutions of Equation (2) 

From Ref. (5), equation (2) can be expressed by a series 
introduced by A. Clebsch. Since * is measured on x-axis, we can 
keep the terms containing only cosn«P . Moreover, there are two 
concentrated forces, hence n is an even number. Then 


u>=/?o(r)+ E P,(r)cOSnr. 

(9) 

in which the coefficients represented by R (r) and R (r) 

o n 

(n=2,4,6, .. .) are the functions of the polar radius "r". (They 

will be abbreviated as R , R ) . When equation (9) is substituted 

o n 

into Equation (2), we obtain the following differential equation 
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( 10 ) 


/ d l 1 d_^\(d^ ?[£. __2!£s_\. 0 , 

(dT' + T > r*A rfr* r f* ^ 


Let us introduce a non-dimensional variable p* r/r b , the we can Z 1 ?— 
solve equation (1) as 


r J£,0‘ + C|lo/9 + r\D a (i l \nf} ( 3 n - 0 ) » 

Rn* r\A n (}* + + rl* l C*&** 1 + (atn>2 ) . 


( 11 ) 

( 12 ) 


When we utilized the boundary conditions, i.e., equations 

(6), (7) and (8), we can obtain the values of the coefficients A ( 

, B , C , D , A , B , C and D . 
o o o o o o o 

For the solution R q , when n * 0, the group of linear 
equations containing the coefficients A Q , B q , C 0 , and D q can be 
derived from the boundary conditions. After the solution is 
obtained, when we assume 

a a * (nD/Prl)A 0f 6„ = (*£)/ P)B tt c t * (nD/ Pr\)C„ d t = (nD/P)D t , 
then "R " can be expressed as 


*. 3 (Pr\/nD){o a + ♦c.ln/? + </ e 0 J ln/0), 

in which all the coefficients can be expressed as 

( o 0 = - b 0 (J\ -Coin/?, -d a f}\\n0it 

c, (1 7 <3+0) 

6c= '2 <1 +#•) ” * 2(1 ♦#«> p 

Cb=-O 0 Pi ^,(l-M)+^i , (l+/ i ) 

( <f 0 = -0.25, 

where p, = r c /r b . 

As to the condition n^i, when we consider 


(13) 


(14) 
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we can obtain the group of linear equations containing the 
coefficients A n , B R , C n , and D n , which can be obtained by solving 
this group of linear equations. Let 

a.«(* D/Pr?")A. t 6.« (*D/Pr;")B„ c.« (*£)/ Pr\')C„ («D/ Pr\)D. . 

the R n can be expressed as 


where 


R » « (Prl/nD)(a,/3 t + bj~ n + c m 0"' + dj"") 


(15) 




n(n + l)(l -*0# 




r»(r» - 1 >(1 


«-(i -«>(»• -i) ♦(«*••■ 


(16) 


Therefore, the axial displacement "w" of any point on the middle 
plane of the flat flexible gear will be 

... _ M r 1 / nfl * *#)./?*• + r ./? 4 41 — 


u> * — — (a, + + c e ln/9 + d t p l \np) + ^ +c«^ 


+ </.^-* 4 *)C0S^]. (17) 

When p= 1 and O', the axial displacement "w" reaches maximum 
"w ", then from equation (17), we can use "w o " to express the 
equation for an axial displacement as 

w * — f(a, + 6 0 /?* + c,ln£ + +djl Mn0) 

(a 0 + 6|)+ C (o. + 6, + c. + </.) L (18) 

♦ 2 <a.0* ♦ b.0- + ♦ d.0’"') ocmfKp] 
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In practical design, w Q for the end harmonic gearing 
transmission flexible gear Is always determined in the design 
specification. Hence, equation (18) is very simple and easy to 
use. 

IV. RADIAL DISPLACEMENT AND TANGENTIAL DISPLACEMENT 
From the first equation of (3), the radial displacement "u" 
can be expressed as 


• ■ - t/(Tr)'*” * j ri ‘l(-wY Jf ’ 


Bui \* 


when we substitute (18) In this equation, and let 

+ £ (o. 6. + c. ♦</.), then after inte- 

gration, we obtain 




»«-<u/J /2)r;'*i*[c/, + 2 E £ U.o <v» 

• ^44.- 


•■ 144 . 

+ £ C U*m<xmn<poo§m<p 


In which 


(19) 


U, . +rf|[(In^)* - fl«# + -i) 

+ M. j + 20c, [i, *2d,lt,0-d,1-cl0-‘, 

a - 

U - ■ (^Tnl 2 " 6 * 0 * + "I* 0 * ♦ <" ♦ 2 )f.c. ♦ 2~/,o.(ln^ - 

+ (fM)[ 2n6c6 ' + "*•*• 4 <" " 2) ♦ 2m/,6.(ln/9 + ^)] 

♦ (aVsjl 2( " + 2)6 * c ' + + 2)rf|C ' + 2( " * 2></,c.(ln^ - j^)] 

B' % ** r 

+ 0T^3>L 2( " ’ 2 )b ' d > + - 2)c/,rf. ♦ 2(» - 2)d.dJ\nB ♦ ^)] 

^•-1 


( 20 ) 


( 21 ) 
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/ 180 




W*Q» I 

(2n - 1) 


0‘*-' 


n'bl „ 
2n+ 1 ^ 


»•- 


+ 


( W 4 2)*C>* «| », | 
(2«Vj) p 


(n- 2)^ 
( 2n - 3) 


0‘f*» 4, gn'Odft.^* 1 


2w(w + 2)a.c < /7 ,.., 
(2r + 1 ) p 


2n(fi-2)fr.</. / Q.,, 4 i 2(n-2)(n + 2 )c„d, „ 
(2n- 1) P 3 P 

- 2n(n- 2)a,d t 0 -2 n(n + 2)6.r ./?» 


i 




nma„ g m ( 

n + m- ] “ 


nna,b m 
n - m - ] p 


/?•*"* ' 

♦ Cw(m + 2)o.c,+ w(m-2)c,o,D . ri- 

_ [n(m - 2)o»d m + m (n + 2)c,6„] , 

n- rn ♦ 1 p 

_ om -*- 1 _ wmf>,6 w - , 

m - n - V ~ n + m + 1 p 

_ C"(”« ♦2)fe ,c „ w (w - 2)d,a m J 
n - ra + ] 

C wQi - 2 )rf.6, ♦ n(m - 2)6. </«] . . 

« 4- m - ] 6 


+ M2)(w 4 2)c, f . , (« + 2)(m - i)c „d m . 

«+m4j ^ n-m+J p 

_ JJL-_2)_<w»2)rf.c . (w-2)(m-2) < /.c. . 

m - n + 3 p n-fm-3 p • 


The integration constant C n (p,f) can be determined by the 
boundary conditions. When P* Pj » snd u * 0, then 

) - ~ u-iri 1 *; 1 [t/,4-2 E 

»•••) •* * 

♦ E U'CM : n<r>+ E E l/^sosarsossirj . 


(’ 2 ) 


(23) 


(24) 
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The tangential displacement v can be determined by the 
second equation of equation (3) t thus 

» * "f[“ *fr(w) “K" + 2rb(“^*) If** 9 (25) 

When we substitute equation (19) and the results of the second 
derivative from the equation (IS) with respect to f, after 
Integration we obtain 


(26) 


where 


V • fCr(/9, <P), 

V. m 2 E ^j- <(/„-(/,.,)« imv + E -<£/.- (/,,)(-—-♦ ~»i o2e*) 

i*ui. " ' * * * 

t-. 1 **"’*-. ,, S ■in(m + n)f . »in(m - n)9? "\ 

*J£ a •* ' * “mSUT > (27) 


in which t) . IJ , and U explain the values of 11 , U and 

con cn cum r on n 

U nm when P-P,. They can be obtained by substituting p*p, 
respectively in the equations (20), (21) and (23). 


V *•- E 

•■I4A- * 

jnn(m»n) g> *| \ 

♦ EE *•”•/?,/?» 2\m-n) ~ 2("» + «) J* 


/ 181 


(28) 


The integration constant C v (8,f) can be obtained by 
considering the following two conditions: when p*p.j, v * 0; and 

when e ■ 0* and 90*, v * 0. The physical meaning of the latter 
is apparent: in order to satisfy the conditi n that the strain in 
the middle plane is zero for any value of "r", when « ■ 0* and 90* 
"v" has to be zero. Therefore, we can Ignore the terms which do 
not satisfy the boundary condition in equations (28) and (27), 
hence , 


Ct( 0 ,V) m **) ■ " 


(29) 
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V. ORIGINAL CURVED SURFACE EQUATION C p 
Similar to studying of the engagement theory for the 
harmonic gearing transmission, we define the middle plane 
equation after the deformation of the flexible gear with the 
initial curved surface equation Based on the displacements 

u, v, and w, we can obtain C p . 

Let us assume represents the equivalent rotating angle of 
the vector generated by the tangential displacement of any 
considered point on the flexible gear, then C p can be expressed 
by the first order approximation as 

( jr«(r»l*i>eoe?,i 

< K ■ ii 

(z-tt r. (30) 


in which, ?>, * 9P ♦ vfr t 0 , 

Equation (30) is a basic equation for studying the end harbor, ic 
gearing transmission geometry. 

VI. CALCULATION RESUlTS AND CONCLUSION 
For our convenience, we can write Equations (18), (19), and 
(26) in the following forms 

! y ■ - n'(urS/r»)p 

i> ■ v'(tvi/r»>i (31 ) 

u; ■ ur'uN* 

In application, for a flexible gear: r c - (0.3 *^,0.4)r^, 

therefore we conducted a series of calculations assuming P^« 0.3 
and 0.4. In the calculations, we select * 0.3, and the terms 
in the series n - 2, 4, 6,... 20. The calculated results u', v', 
and w', from a FfeLIX C-256 computer, are shown in Figure 2 (P^ ■ 
0.3) and Figure 3 ( 0.4), in which the solid lines represent 
the values of w'; dotted lines, v'; and dot-dash lines, u'. From 
Figure 2 we can see whatever p varies, the maximum values of w' 
and u' always occur at t- 0*, yet the values of "v" 1 always vary 
with p. For example, when p » 1.0, the maximum v' occurs at 
around t* 15*; and P* 0.6, 25*. Figure 4 shows the curves of 
the maximum values if u', v', and w' vs p. 
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When pal, the values of u, v, and w are critical for the 
engaging analysis of the end harmonic gearing transmission. 
Hence, for the designer's applications, we calculated the values 
(Table I) of u', v', and w' under the conditions of p^= 0.3 and 
0.4 respectively when p= 1. 




Figu? 3 Values of u', v' and w' , when p-= 0.4. 





Figure 4 p-j * 0.3, u',v' and w' vs p. 



Figure 5 Theoretical Analysis of Actual Measurement 
( P 1 = 0.37, p = 1) 

In order to verify the accuracy of this calculation method, 
we conducted the measurement of a flat flexible gear (with the 
following dimensions: number of teeth z = 268, gear module m = 
0.8mm, tooth width b = 13mm, thickness v * 2.0mm, r = 40mm) 

V 

under the condition that w Q = 1.6mm. Figure 5 shows the results. 
In this figure, the solid line represents the calculated values, 
and the circles represent the measured data (each datum point is 
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the average of five measurements). The measurement shows calculated 
value and the measured value exceed no more than 0.015mm. 

Hence, the error is within the tolerance and the two values are 
very close. It is concluded that the calculating method 
recommended in this paper can meet the precision requirement of / 1 83 
the end harmonic gearing transmission analysis, and is a feasible 
method. 

Table I. Values of u', v', and w' , when p= 1.0 


*>• 



fix = 0.3 

A = 0.4 


M * 

v ' 

w' 

! ✓ 

v ’ 

w' 

0* 

0.8321728 

0.0000000 

1.0000000 

0.9775204 

0.0000000 

1.0000000 

6* 

0.7965759 

0.0764656 

0.9798877 

0.9260779 

0.0868258 

0.9749692 

10* 

0.7017223 

0.1314390 

0.9235046 

0.7916194 

0.1416512 

0.9054893 

15* 

0.5757025 

0.1561112 

0.8410518 

0.6198189 

0.1541650 

0.8059808 

20* 

0.4468286 

0.1556362 

0.7449075 

0.4541346 

0.1364796 

0.6937554 

25* 

0.3325486 

0.1406637 

0.6453400 

0.3179190 

0.1071461 

0.5827968 

30* 

0.2385098 

0.1192367 

0.5481229 

0.2150134 

0.0777533 

0.4803274 

35’ 

0.1636258 

0.0947998 

0.4648851 

0.1396382 

0.0513040 

0.3873771 

40* 

0.1052356 

0.0684704 

0.3653760 

0.0849351 

0.0271207 

0.3021504 

46* 

0.0615529 

0.0417018 

0.2798654 

0.0484686 

0.0051485 

0.2235869 

60’ 

0.0314360 

0.0172052 

0.2001617 

0.0217469 

-0.0128444 

0.1529013 

66* 

0.0131310 

-0.0020920 

0.1269178 

0.0080536 

-0.0247142 

0.0926124 

80’ 

0.0038425 

•0.0146868 

0.6081082 

0.0019453 

-0.0299428 

0.0443598 

66* 

0.0004946 

-0.0209949 

0.0180181 

0.0001274 

-0.0299427 

0.0074478 

TO* 

0.0006224 

-0.0224486 

-0.2250659 

0.0004063 

-0.0267263 

-0.0207643 

76* 

0.0025764 

-0.0203651 

-0.0546138 

0.0018714 

-0.0217060 

-0.0430344 

80* 

0.0051916 

-0.0156142 

-0.0792365 

0.0033314 

-0.0154536 

-0.0604660 

85* 

0.0073893 

-0.0084348 

-0.0946413 

0.0047949 

-0.0081088 

-0.0721371 

90* 

0.0082656 

- 0.0000000 

-0.1002553 

0.0053911 

0.0000000 

-0.0763083 
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